Maximum black-hole spin from quasi-circular binary mergers 
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Black holes of mass M must have a spin angular momentum S below the Kerr hmit {% = S/M'^ < 
1), but whether astrophysical black holes can attain this limiting spin depends on their accretion 
history. Gas accretion from a thin disk limits the black-hole spin to Xgas ^ 0.9980 ± 0.0002, as elec- 
tromagnetic radiation from this disk with retrograde angular momentum is preferentially absorbed 
by the black hole. Extrapolation of numerical-relativity simulations of equal-mass binary black-hole 
mergers to maximum initial spins suggests these mergers yield a maximum spin Xeg lS 0.95. Here we 
show that for smaller mass ratios q = m/M <^ 1, the superradiant extraction of angular momentum 
from the larger black hole imposes a fundamental limit xiim ^ 0.9979 ±0.0001 on the final black-hole 
spin even in the test-particle limit (g — > 0) of binary black-hole mergers. The nearly equal values of 
Xgas and xiim imply that measurement of supermassive black-hole spins cannot distinguish a black 
hole built by gas accretion from one assembled by the gravitational inspiral of a disk of compact 
stellar remnants. We also show how superradiant scattering alters the mass and spin predicted by 
models derived from extrapolating test-particle mergers to finite mass ratios. 



I. INTRODUCTION 

Supermassive black holes (SBHs) reside in the centers 
of most large galaxies. While a iew nearby SBHs can be 
detected by their gravitational influence on surrounding 
stars, the majority of SBHs are observed electromagnet- 
ically as active galactic nuclei (AGN). The same accre- 
tion flows that release energy to power AGN also supply 
energy and angular momentum to the black holes them- 
selves increasing their mass M and spin 5' [43]. By defi- 
nition, black holes possess an event horizon from within 
which nothing can escape to future null infinity. Black 
holes described by the Kerr metric [J only possess an 
event horizon for x = S/M"^ < 1, setting a fundamental 
upper limit on a black hole's possible spin. Whether or 
not astrophysical SBHs saturate or even exceed this Kerr 
limit is an important test of general relativity. 

Black-hole spins also probe their assembly history. 
SBHs grow both by gas accretion and mergers driven by 
the gravitational inspiral of binary companions. These 
two growth mechanisms may supply mass and angular 
momentum to the SBHs in different ratios, allowing mea- 
surements of black- hole spin to distinguish between them. 
Reverberation mapping of the iron Ka line in AGN x-ray 
spectra has been proposed as just such a means of mea- 
suring SBH spins [2]- XMM-Newton observations ana- 
lyzed with this technique have been used to constrain 
the spin of the SBH hosted by the Seyfert 1.2 galaxy 
MCG-06-30-15 to X = 0.989i|]™^ at 90% confidence [3]. 
A closer examination of how SBHs acquire their spins is 
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thus of both theoretical and observational importance. 

According to Bardeen [4] , a nonspinning black hole can 
attain the maximum Kerr spin x = 1 after accreting a fi- 
nite mass of test particles freely falling from its innermost 
stable circular orbit (ISCO). However, material on circu- 
lar orbits with radii greater than that of the ISCO can- 
not be accreted unless it has some mechanism to shed its 
excess angular momentum. In a standard geometrically 
thin, optically thick accretion disk "5^, this mechanism is 
viscous stress within the disk that also heats it locally 
and transports energy outwards. This heating produces 
an energy fiux at the disk's surface that will be radiated 
in all directions. A small fraction of these radiated pho- 
tons do not escape to infinity but are instead captured 
by the black hole itself. Photons with negative angular 
momentum with respect to the black hole's spin have a 
larger capture cross-section than those with positive an- 
gular momentum [6], implying that the accreted radia- 
tion will counteract the directly advected material which 
always acts to spin the hole up for x < 1- These two 
sources of angular momentum cancel for black-hole spins 
Xgas — 0.998, with black holes of greater spins spinning 
down to this value after accreting a mass m ~ 0.05M [7]. 

Different accretion flows will supply energy and an- 
gular momentum to the black hole in different ratios, 
altering the value of the limiting spin Xgas- Advection- 
dominated accretion flows (ADAFs) do not cool effi- 
ciently, and therefore a fraction / > of the gravita- 
tional energy dissipated prior to the ISCO is advected 
by the black hole rather than radiated to infinity. Vis- 
cous stresses in these flows proportional to the Shakura- 
Sunyaev parameter a [5 also reduce the speciflc angular 
momentum of the accreted material below its value in 
the vacuum at the ISCO. The former effect increases M 
in the denominator in the definition x = S/M"^ of the di- 



mensionless spin, while the latter effect reduces S in the 
numerator [8j . The magnetorotational instability [9] also 
torques gas at the ISCO, and in addition can launch jets 
which further limit the spin as shown in magnetohydro- 
dynamic simulations [TU] . Analytic fits to the ADAF sim- 
ulations of Popham and Gammie [S\ suggest that black 
holes can spin up to xadaf — 0.96, and that the in- 
clusion of jets calibrated by the magnetohydrodynamic 
simulations of Gammie et al. (TUj reduces this limit to 
Xjets — 0.93 [11]. We see that limits on black-hole spin 
depend greatly on the nature of the accretion flow, and 
that no model-independent constraints exclude the limit 
Xgas — 0.998 for thin disks first obtained by Thorne [7]. 

The maximum spin of black holes produced in binary 
mergers is also greatly uncertain. The most accurate 
numerical-relativity simulations have been of equal-mass 
q = milrrix = 1 black holes with spins xi — X2 aligned or 
antialigned with the orbital angular momentum. These 
simulations found x = 0.68646 ±0.00004 for nonspinning 
binaries [l2], and x = 0.547812 ± 0.000009 for equal- 
mass binaries with initial spins Xi = X2 = —0.43757 |13j . 
The spins of such aligned binaries do not precess prior to 
merger, removing one complication in determining their 
magnitude and direction. We will therefore restrict our 
analysis to aligned binaries in this paper, though with 
modest additional work it could be extended to binaries 
of arbitrarily oriented spins. 

Although these numerically determined spins are re- 
markably precise, two significant obstacles still prevent 
the determination of final black-hole spins from more 
generic mergers. The first is the inability of the com- 
monly used conformally flat Bowen-York initial data to 
adequately approximate initial binary black holes with 
spins greater than Xi — 0.93 [H]. An extrapolation of 
simulations with Xi ^ 0.9 suggests that the merger of 
equal-mass black holes with maximal initial spins aligned 
with their orbital angular momentum yields a final black 
hole with spin x = 0.951 ± 0.004 [15]. However, it re- 
mains to be proven that such a smooth extrapolation in 
initial spin holds all the way to the Kerr limit. 

The second obstacle to determining spins from generic 
mergers is the increased computational resources needed 
to simulate mergers with mass ratios q < 1- This in- 
creased computational demand results from the need for 
more closely spaced numerical grid points to resolve the 
horizon of the smaller black hole, shorter timesteps ow- 
ing to the shorter light propagation time between grid 
points, and longer simulations to capture the larger num- 
ber of orbits per unit increase in orbital frequency prior 
to merger. The smallest mass ratio that has currently 
been simulated is q — 0.01 [IB]; only a handful of sim- 
ulated mergers with g < 0.1, all of initially nonspinning 
black holes, have been published [T3{19]. New numer- 
ical techniques will be required to make much progress 
beyond this point. 

Relativists have attempted to surmount these obsta- 
cles by inventing fitting formulas that are functions of 
g, xii ^-nd X2, calibrating the coefficients in these for- 



mulas with existing simulations, and then extrapolating 
them to higher spins and lower mass ratios than can cur- 
rently be simulated. This approach is very effective in the 
region of the parameter space {g,XijX2} near the simu- 
lations with which the fitting formulas were calibrated, 
but can break down outside this region. The probable 
reason for this is that all three parameters listed above 
equal unity for equal-mass, maximally spinning, aligned 
mergers, and thus polynomials in these parameters will 
converge slowly if at all. Even the symmetric mass ratio 
T] = mim2/(mi + 1712)^ = q/{l + q)^ gets as large as 1/4 
for equal-mass mergers, suggesting that polynomials in 
this parameter will not converge quickly either. At the 
time this paper was written, simulations were restricted 
to the region 0.1 < q < 1 (0.0826 < tj < 0.25) which pro- 
vided a short lever arm over which to calibrate terms with 
different q dependence. The degeneracy between terms 
makes predicting final spins in the test-particle {q — >■ 0) 
limit very uncertain. For example, several fitting formu- 
las predict final spins above the Kerr limit for maximally 
spinning aligned mergers with mass ratios as modest as 
q < 0.25. The publication of a simulation with q = 0.01 
|16) during the preparation of this paper shows the rapid 
progress towards numerical simulations of extreme-mass- 
ratio mergers, but much work remains before such simu- 
lations are available for generically spinning black holes. 

Fortunately our analytical knowledge of the geodesies 
of the Kerr metric can help us understand binary black- 
hole mergers in the test-particle limit. These mergers 
can be broken down into three stages: inspiral, plunge, 
and ringdown. During the inspiral stage, the orbit of 
the test particle adiabatically evolves through a series 
of geodesies with successively lower energy and angular 
momentum as these quantities are radiated away through 
the emission of gravitational waves. This gravitational 
radiation circularizes initially eccentric orbits while the 
evolution is post-Newtonian [5D] , implying that for many 
though possibly not all astrophysical mergers the orbit 
will have circularized before the test particle reaches the 
ISCO. After this point, during the brief plunge stage, 
the test particle rapidly falls into the event horizon of 
the larger black hole. Finally, the quasinormal modes 
excited during the merger ring down as the newly formed 
black hole settles into its final Kerr configuration. In the 
test-particle limit, the energy and angular momentum 
radiated during the inspiral stage scales linearly with q, 
while that radiated during plunge and ringdown scales as 
a higher power in the small parameter q. This suggests 
that one can reasonably predict the mass and spin of the 
final black hole by equating these quantities to the energy 
and total angular momentum of the binary at the ISCO. 

This prediction assumes that during the inspiral all of 
the gravitational waves are radiated outwards to infinity 
rather than downwards to the larger black hole's event 
horizon. For highly spinning black holes, these downward 
gravitational waves will be supcrradiantly scattered, ex- 
tracting energy and angular momentum from the larger 
black hole. Individual modes can be amplified as much as 



138% by a maximally spinning black hole [3T]; the total 
energy flux radiated to infinity by a test particle at the 
ISCO will be amplified by 12.9% for a black hole of spin 
X = 0.999 [25]. Although this amplification decreases 
rapidly as a function of the orbital radius, the total en- 
ergy and angular momentum extracted throughout the 
inspiral still scale linearly in q and will therefore remain 
the dominant factor affecting the final mass and spin af- 
ter the energy and angular momentum advected with the 
test particle itself during merger. 

The primary goal of this paper is to determine how 
much the superradiant scattering of gravitational waves 
during the inspiral reduces the maximum spin achievable 
by binary black-hole mergers below the Kerr limit first 
predicted by Bardeen [1]. In Sec. In] we will review the 
previous analytic predictions for the final mass and spin 
that serve as our starting point. Our method for calcu- 
lating the energy and angular momentum extracted from 
the spinning black hole will be described in Sec. |III[ T he 
results of this calculation are presented in Sec. |IV| In 
Sec. Ivlwe compare our model to others in the literature, 
extrapolate it to comparable-mass mergers and compare 
with numerical-relativity simulations, and suggest how it 
might be combined with fitting formulas to produce an 
approximation valid over the entire domain < q < I. A 
brief summary and a few final remarks on the astrophys- 
ical implications of our analysis are offered in Sec. |VI[ 



II. TEST-PARTICLE MERGERS 

Hughes and Blandford f53] (hereafter HB) recognized 
that in the test-particle limit, the energy and angular mo- 
mentum radiated during plunge and ringdown was much 
less than that released during the inspiral stage. They 
predicted the mass Mf and spin Sf of the final black hole 
would therefore be given by 



5'/,HB = rnim2LiscoiXi) 



mfxi 



(la) 
(lb) 



Here, -Eiscolx) is the energy per unit mass of a test par- 
ticle on the equatorial ISCO of a Kerr black hole with 
dimensionless spin x- ^isco(x) is the corresponding di- 
mensionless orbital angular momentum. The dimension- 
less final spin is simply Xf — Sf/M'i. This formula is 
exact in the test-particle limit, and for g — > 0, xi — )• 1 
correctly reproduces the Bardeen result 



9Xf, 



HB 



dq 



^isco(l)-2£;isco(l)=0 



(2) 



Maximally spinning black holes cannot be spun up above 
the Kerr limit. 

Although Eqs. (II]) are exact in the test-particle limit, 
they are manifestly asymmetric in the black-hole labels 
"1" and "2" . For example, they ignore altogether the spin 
X2 of the smaller black hole. While HB boldly extrap- 
olate their formula to g = 0.5, Buonanno, Kidder, and 



Lehner [M] (hereafter BKL) realized that a symmetrized 
version would more accurately describe comparable-mass 
mergers. They proposed 

-^/,BKL = TOi + m2 , (3a) 

•S/^BKL = "ll?™2ilSCo(X/) + TOiXl +"12X2 ■ (3b) 

In addition to the obvious improvement of including the 
second spin X2, BKL made the inspired choice of using 
the dimensionless spin x/ of the final black hole to calcu- 
late the orbital angular momentum of the initial binary at 
the ISCO. Although counterintuitive at first, this choice 
correctly captures the total energy and angular momen- 
tum of the system which are assumed to be conserved 
after the binary reaches the ISCO. 

Equations. (Jsj) are far more successful at predicting the 
final spin from equal-mass mergers than they have any 
right to be. They predict that equal-mass nonspinning 
black holes will merge to yield a final black hole with spin 
X/,BKL — 0.663 quite close to the numerically determined 
value XNR = 0.68646 ± 0.00004 |I2]. Equal-mass, max- 
imally spinning black holes are predicted to produce a 
final spin x/,bkl = 0.959 which is also surprisingly close 
the numerically extrapolated value xnr = 0.951 ± 0.004 
[T5] . This latter agreement however is an artifact of over- 
estimating the final mass in Eq. ( 3a ) by failing to account 
for the energy radiated during the inspiral. Overesti- 
mating the denominator in the expression x/ — Sf/Alj 
leads to an underestimate of the final spin. In the limit 
q — >■ 0, xi ~^ 1; tttis underestimation implies 



<9X/,BKL 

dq 



iisco(l)-2 = 2(3-i/2-l)<0. (4) 



Maximally spinning black holes are artificially found to 
be spun down by test-particle mergers. In this model 
black holes can only be spun up by test particles to the 
fictitious limit xiim = 0.948 at which dxf/dq = 0. 

Kesden [5S] (hereafter K) sought to remedy this by 
replacing Eq. (3a) with 



Mf^K = M - fi[l - EiscoiXf)] 



(5) 



where M = mi + m2 is the sum of the initial masses, 
fi = mim2/M is the reduced mass of the binary, and the 
energy per unit mass £'isco(x/) is evaluated using the fi- 
nal dimensionless spin in the spirit of BKL. This formula 
retains the desired symmetry of Eqs. ([3J) under exchange 
of black-hole labels, but also reduces to Eqs. ([ij) in the 
test-particle limit thereby preserving the Bardeen result 
that black holes can be spun all the way up to the Kerr 
limit by test-particle mergers. It predicts that equal- 
mass, nonspinning black holes will merge into a hole 
with final spin x/.k,ns = 0.687 in near miraculous (and 
probably coincidental) agreement with numerical simula- 
tions, but also predicts an uncomfortably large final spin 
X/,K,s — 0.9988 for maximally spinning aligned mergers. 
Equations (3b) and (l5|) are not the unique choice that 
possesses the desired symmetry and limiting behavior. 



For example, replacing Eq. (3b) with 

Sf,K = f^MfLisco{Xf) + mlxi+mlx2 (6) 

also maintains these properties, but predicts different re- 
sults X/,K,NS — 0.675 and X/,k,s = 0.9909 when extrapo- 
lated to equal-mass mergers. For our purpose, in this 
paper of calculating the maximum final spin for test- 
particle mergers, either Eqs. (la) or Eq. ([5| paired with 
either Eq. ( 3b ) or Eq. ^ can serve as a suitable starting 
point. 

III. SUPERRADIANT SCATTERING 

The gravitational radiation emitted by an inspiraling 
test particle is fully described by the complex Weyl scalar 



Since r is negative as are the fluxes (dE/dt) 



rad 



{dL,/dt) 



V'4 



-C aP'ysn"' m rC ra 



(7) 



where Cap-yS is the Weyl curvature tensor and n^ and 
ifi^ are elements of a Newman-Penrose tetrad of null 4- 
vectors [55]. Teukolsky [57] decomposed tp4^ into multi- 
pole moments 



■04 = (r— mcos( 



duj ^ Rw. 



(r)_25r-(0)e™*e 



(8) 
and showed that the evolution of the radial modes 
RimLoij") is governed by an ordinary second-order differ- 
ential equation sourced by the test particle's contribution 
to the stress-energy tensor. Sasaki and Nakamura j28.' de- 
rived a linear transformation of Rimto {i") that greatly fa- 
cilitates the solution of this differential equation. To solve 
this equation ourselves, we used the GREMLIN (Gravita- 
tional Radiation in the Extreme Mass ratio LIMit) code 
written and generously provided to us by Scott Hughes. 
We adopt his notation and closely follow his treatment 
presented in |29j to calculate the energy and angular mo- 
mentum extracted from the black hole by superradiant 
scattering during the inspiral. 

The Teukolsky-Sasaki-Nakamura formalism allows one 
to calculate the rates {dE/dt)^!^^^ and {dLjdt)l''j^^^ at 
which energy and angular momentum are radiated down 
the event horizon at r+. The total energy i^sR £^nd an- 
gular momentum JgR extracted from the black hole as 
the test particle inspirals from infinity to r are integrals 
of these fluxes 



^sr(x,'") 
Jsviix^r) 



where the radial velocity f during an adiabatic, quasicir- 
cular inspiral is given by 



dEY'' dr' 


(9a) 


dLA'""^ dr' 
dt /^_,^ r 


(9b) 



r{x,r) = 



dEV"^ UE\~^ 



dt J^^^ \ dr 



(10) 



rad 



and 

for high spins, £'sr and Jsr are defined 
to be positive when energy and angular momentum are 
extracted from the black hole and negative when they 
are absorbed. Adding the contributions of Eqs. ^ to 
the right-hand sides of Equations ([s]) and ( jSb) ) yields 



Mf = 
Sf = 



Xf 



M - ^[1 - £^isco(x/) + £'sr(x/, risco)] (Ha) 
TOiTO2[iisco(x/) - JsKiXf^nsco)] (lib) 



-mlxi 



»7l2X2 



SflM] 



(lie) 



for our revised prediction for the final black-hole mass 
and spin. With a slightly different extrapolation of Sf to 
large mass ratios, we can modify Eq. (rol to obtain 

S'f = ^iMf[Llsco{Xf) - JsK{Xf,nsco)] (12) 
+mlxi + mlx2 ■ 



In the test-particle limit q ^ 0, Eq. (11a) combined with 
either Eq. ([Tib]) or Eq. ([l2| gives 



-g^iXi,r) = —^ (Xi)--^(xi,0 (13a) 



t^Xisco 
dq 

%SR 



(Xi) -^ iisco(xi) -2xi£^isco(xi) (13b) 
{Xi,r) -^ JsR{xi,r)-'2xiEsn{xi,r) .{13c) 



At the limiting spin xi = Xi™, dxfixi,nsco)/dq = 0, 
implying that black holes cannot be spun up beyond xum 
by test-particle mergers. 



IV. RESULTS 

In the top panel of Fig. [T] we show the dimensionless 
energy £'sr(Xj^) and angular momentum Jsr(X:^) ex- 
tracted from a black hole of spin xiim = 0-9979 by super- 
radiant scattering as a test particle inspirals from infinity 
to a radius r. According to Eq. (13c), these two quan- 
tities determine the spin decrease dxsn/^q due to the 
superradiant scattering of downward-going gravitational 
waves. This spin decrease is shown for spins x = 0.996, 
0.997, xiim, and 0.999 in the bottom panel of Fig.[l] The 
increase in spin dxisco / 9q when the test particle itself is 
accreted by black holes of the same spins is shown by dot- 
ted horizontal lines. At the limiting spin xiim, the total 
spin decrease dxsR/dq evaluated at the ISCO precisely 
cancels the spin increase dxisco/dq leaving the dimen- 
sionless spin X unchanged. This is seen explicitly by the 
intersection of the solid black curve and dotted black line 
at nsco in the bottom panel of Fig. [l] The intersection 
of dxsn/dq and dxisco/dq at xi = Xiim, r = nsco 
implies that the total spin change dxf/dq vanishes by 
Eq. (13a). Note that both the spin angular momentum 
Sf and mass Mf do increase as a result of the merger. 



but in just the right ratio as to preserve Xf = Sf /Mi 
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FIG. 1: Top panel: The dimensionless energy Esn (thin hne) 
and angular momentum Jsr (thick hne) extracted from a 
black hole with the limiting spin xiim ~ 0.9979 as the test 
particle inspirals from infinity to radius r. Bottom panel: 
The change in final spin per unit test-particle mass dx/dq. 
The dotted horizontal lines show the spin increase dxisco/dq 
when the test particle falls from the ISCO into a black hole 
of spin X — 0.996 (green), 0.997 (blue), xiim (black), and 
0.999 (red). The curves Oxsn/dq show the spin decrease as 
superradiant scattering extracts energy from black holes with 
spins of the corresponding color. Only at x = Xiim do the 
lines and curves intersect at risco, indicating that the spin 
remains unchanged by the merger. 



The spin increase dxisco/dq is a monotonically de- 
creasing function of xi, while the spin decrease dxsn/dq 
evaluated at nsco monotonically increases with xi ■ This 
implies that the total spin change dxf/dq monotonically 
decreases with xi, intersecting zero at xi = Xiim- Black 
holes with xi < Xiim are spun up by mergers with test 
particles on circular equatorial orbits, while black holes 
with spins greater than xiim are conversely spun down. 
Black holes are never spun down [dxf{xij''"iSCo)/9q > 
0, Vxi] when superradiant scattering is neglected, repro- 
ducing the original Bardeen |4] result that the Kerr limit 
could be saturated. 

Only one curve for i?sR and one for Jisco were pre- 
sented in the top panel of Fig. [l] as the difference in 
these quantities as functions of the initial spin could not 
be distinguished in the logarithmic plot needed to depict 
the five orders-of magnitude of their evolution. In Fig. [2] 
we zoom in on the region near the ISCO to display this 
spin dependence. The total energy and angular momen- 
tum extracted are nearly independent of spin all the way 



FIG. 2: The same quantities Esn, Jsr, and dx/dq presented 
in Fig.fllas functions of the orbital radius r. We zoom into the 
region near the ISCO to better display how these quantities 
evolve near merger. The colors and styles of the curves are 
the same as those in Fig. [T] 



down to r = 2mi, deep within the relativistic region. 
Only very close to the ISCO do the curves diverge, with 
the higher values of Esr and Jsr for more highly spin- 
ning black holes coming in nearly equal parts from the 
higher fluxes at fixed radii and the decreasing value of 
nsco in the lower bound of the integrals in Eqs. ([9]). 

Using this refined model for test-particle mergers, we 
can recalculate how the dimensionless black-hole spin x 
and mass M evolve after the black hole has accreted a fi- 
nite mass AM of test particles gravitationally inspiraling 
inwards from large radii. Following Thorne |7], we find 
that this evolution is governed by the coupled first-order 
differential equations 



dx 
dAM 

dM 
dAM 

dJ 
dAM 



1 dJ 



2x dM 



M2 dAM M dAM 
Eiscoix) - EsB.{x,nsco{x)) 



(14a) 
(14b) 



M[Liscoix) ~ JsR{x,nsco{x))] ■ (14c) 



In Fig. I3J we see how the spin x evolves for different val- 
ues of the initial spin xo (where the curves intersect the 
y axis AM/Mq — 0). The black curves show the predic- 
tions of our model, while the blue curves show those of 
Bardeen [4] without the scattering of gravitational waves 
by the black hole. As expected, in our model black holes 
with initial spins x ^ Xum can only be spun up to this 
limiting value no matter how much mass they accrete. 
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FIG. 3: The spin x of 8- black hole initially of mass Mo after 
binary mergers with a mass AM of test particles on quasi- 
circular equatorial orbits. The black (solid) curves show the 
predictions of our model for different initial spins, while the 
blue (dashed) curves show the predicted spins in the absence 
of scattering. Spins in our model asymptote to xiim — 0.9979 
as shown by the horizontal dotted line, while without scatter- 
ing black holes can spin all the way up to the Kerr limit. 



More highly spinning black holes are spun down to near 
this limit after accreting > 10% of their initial mass. 
By contrast, in the Bardeen model black holes reach the 
Kerr limit x = 1 after accreting a jinite mass AAf of test 
particles. For initially nonspinning holes, this mass is 



AM = 3Mo(sin" 



-il/3). (15) 



Radiatively efficient black holes have a luminosity 



dAM 



e- 



dt 



(16) 



where e — 1 — Eiscq{x) is the radiative efficiency. In 
Fig. J4| we show the efficiency e of the black holes in 
Fig. |3[ as a function of the total mass AM of test parti- 
cles they have accreted. If the superradiant scattering of 
gravitational waves is neglected, the black holes spin up 
to the Kerr limit x = 1 after which they are capable of 
converting a fraction e = 1 — l/\/3 = 0.423 of accreted 
mass into radiant energy. By contrast, accounting for 
superradiant scattering reduces the maximum black-hole 
spin to xiini = 0.998. Although this limiting spin is only 
0.2% below the Kerr limit, the corresponding limiting ra- 
diative efficiency enm = 1 — ^isco(xiim) = 0.320 is 24% 
below that of a maximally spinning black hole. The in- 
ability of astrophysical black to radiate more efhciently 



FIG. 4: The radiative efficiency e — 1 — £isco(x) of the 
accreting black holes shown in Fig. [3] In the Bardeen model 
(dashed blue curves), the black holes spin up to the Kerr 
limit X = 1 after accreting a finite mass A A/. The radiative 
efficiency of a maximally spinning black hole is e = 1 — 1/\/3 = 
0.423. In our model (solid black curves) which includes the 
superradiant scattering of gravitational waves, the black hole 
spins asymptotically approach Xiim- A black hole with this 
spin has the substantially lower efficiency eum = 0.320. 



than Xiim will affect efforts to estimate black-hole masses 
from their luminosities, such as that by Soltan |30j to 
constrain the mass function of supermassive black holes 
from the quasar luminosity function. 



COMPARISON WITH NUMERICAL 
RELATIVITY 



How do the predictions of our model compare with 
those of other published models? In Fig. [5] we again 
show dxf /dq, this time as a function of the initial spin 
Xi- The value of Xi for which a given curve crosses the 
horizontal dotted line dxf/dq = determines the max- 
imum spin Xiim to which a black hole can be spun up 
by test-particle mergers. We see that the superradiant 
scattering of gravitational waves produced during the in- 
spiral reduces xiim from the Kerr limit as predicted by 
Bardeen [3], HB [53], and K [25 and shown by the blue 
(short-dashed) curve to Xi — 0.9979 shown by the black 
(solid) curve. The red (long-dashed) curve shows the pre- 
dictions of the BKL [25 model described in Sec. In] As 
discussed previously, this model artificially reduces xiim 
to Xi = 0.948 by neglecting the spin dependence of the 
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FIG. 5: The change in final spin per unit test-particle mass 
dxf/dq as a function of the initial spin Xi- The solid black 
curve shows the predictions of this paper, with spin down 
(dXf /dq < 0, horizontal dotted line) possible for spins Xi ^ 
Xiim (vertical dotted line). The short-dashed blue curve shows 
how this result changes when superradiant scattering is ne- 
glected as in the models of Bardeen [4] , HB [23] , and K [25] . 
The long-dashed red curve gives the BKL [21] test-particle 
prediction, while the short dash-dotted green and long dash- 
dotted magenta curves are the predictions of the AEI jSlj and 
FAU [321 fitting formulas. 



final mass. The green (short dash-dotted) and magenta 
(long dash-dotted) curves are the predictions of the AEI 
[3T] and FAU [32 "fitting formulas" for the final spins. 

The fitting-formula approach proposes a specific func- 
tional form for the dependence of the given quantity on 
the parameters {?, XI1X2} and a small number of con- 
stant coefficients. These coefficients are adjusted until 
the fitting formula best reproduces the results of a sam- 
ple of numerical simulations. Once the fitting formula has 
been calibrated in this manner, it should be able to pre- 
dict the results of any future simulations. The AEI and 
FAU fitting formulas were calibrated with an extensive 
sample of numerical simulations of binary mergers with 
varying mass ratios and spins, none of which produced a 
black hole whose final spin exceeded the Kerr limit. Yet 
Fig. [5] shows that both curves have dxf/dq > even for 
Xi = 1, implying that test-particle mergers can spin black 
holes above the Kerr limit. One possible explanation of 
this unphysical prediction is that the fitted coefficients 
were misdetermined because of the limited range of mass 
ratios in the sample of simulations with which they were 
calibrated. A second explanation is that the proposed 





Initial spins Xi 


Xi 


-0.43757 ll3J 


[12] 


1 [13 




Final spins Xf 


NR 


0.547812 ± 0.000009 


0.68646 ± 0.00004 


0.951 ± 0.004 


KLP 


0.520861 


0.686354 


0.996439 


KLP' 


0.509269 


0.674197 


0.986947 


K [25] 


0.521153 


0.687036 


0.998805 


BKL [H 


0.505148 


0.663086 


0.959107 


AEI [31] 


0.546646 


0.686460 


0.961491 


FAU [32] 


0.548602 


0.6860 


0.9540 


BK [33] 


0.547562 


0.6893 


0.9504 




Final masses Mf/M 


NR 


0.961109 ±0.000003 


0.95162 ±0.00002 




KLP 


0.979028 


0.974530 


0.920918 


KLP' 


0.979268 


0.974984 


0.932995 


K [25] 


0.979039 


0.974565 


0.916181 


FAU M 


0.962877 


0.9515 


0.9255 


BK [33] 


0.964034 


0.9530 


0.9009 



TABLE I: A comparison between the final spins and masses 
determined by numerical simulations and those predicted by 
various fitting formulas. All three simulations begin with 
equal-mass binaries whose spins are aligned or antialigned 
with the orbital angular momentum and have magnitudes 
given in the first row. These simulations are described more 
fully in the references provided. The next 8 rows give the 
final spin Xf ^s determined by NR or predicted in the refer- 
enced papers. KLP are the predictions of this paper. The 
final 6 rows give the corresponding predictions for the final 
mass Mf/M. 



functional forms are inadequate to describe the mass- 
ratio or spin dependence of the final spin for any choice of 
coefficients. Further numerical simulations, particularly 
at smaller mass ratios, are needed to determine which of 
these explanations is correct. This problem is not just 
restricted to test-particle mergers; the AEI and FAU for- 
mulas predict final spins above the Kerr limit for mass 
ratios as large a,s q ~ 0.283 and q — 0.2434, respectively. 

Although the AEI and FAU fitting formulas break 
down in the test-particle limit, by design they agree 
closely with simulations of comparable-mass mergers like 
those with which they were calibrated. In Table |l] we 
compare our model, extrapolated to equal-mass mergers 
according to Eqs. (11), with both the fitting formulas 



and the most accurate numerical simulations of equal- 
mass {q = 1), equal-spin (xi = X2 = xOj aligned merg- 
ers. Both the nonspinning (xi = 0) [12] and antialigned 
(Xi = —0.43757) [13] simulations were performed with 
the spectral-methods code developed by the Caltech- 
Cornell group. The Xi = 1 numerical result [I^ listed 
in the third column is an extrapolation to maximal spins 
of a series of simulations produced with the BAM finite- 
difference code used by the Jena and FAU groups. 

The first row gives the calculated values of the final 
spin Xf for these three numerical-relativity (NR) simu- 



lations. The next 7 rows give the final spins predicted 
for these binaries by many different published formu- 
las. The second row lists the predictions of this paper, 
Kesden-Lockhart-Phinney (KLP) . Our prediction for the 
nonspinning merger agrees with this numerical result to 
about 10^"'. An agreement this good between a test- 
particle extrapolation and an equal-mass simulation can 
only be a coincidence, as is demonstrated by the third row 
KLP'. Here we have substituted S'f from Eq. ([l2| for Sf 



from Eq. (lib I in our prediction Xf — Sf /Mi. Although 



S'f has the same symmetries and limiting behavior as Sf , 
the prediction changes by about 1% when extrapolating 
all the way to equal masses. Agreement beyond this accu- 
racy must be considered coincidental unless we discover 
a fundamental reason to prefer Sf over S'f. 

The fourth row lists the predictions of Kesden [5S], 
which was the starting point for this paper but neglected 
the superradiant scattering of downward-going radiation. 
Comparing the rows labeled KLP and K, we see that this 
scattering only reduces the final spin by about 0.25% even 
for binaries initially spinning at the Kerr limit. Although 
this effect seems negligible, it is potentially detectable 
since several gravitational-wave observables depend very 
sensitively on x near the Kerr limit. We shall elaborate 



on this in the final discussion in Sec VI Rows 5 through 8 
list the predictions of other published formulas. We have 
already discussed the BKL model adequately; its close 
agreement with the extrapolation to maximally spinning 
binaries is a fortuitous coincidence. The AEI, FAU, and 
BK [33 fitting formulas all do excellent jobs of reproduc- 
ing the numerical results, though the exact agreement of 
the AEI formula with the nonspinning simulation results 
from this simulation being included in the set with which 
this formula was calibrated. 

The final 6 rows of Table |l] show how our predicted fi- 
nal masses Af^, extrapolated from the test-particle limit 
as before, compare with the numerical simulations and 
numerically calibrated fitting formulas. Our predictions 
overestimate Mf compared to the simulations and fitting- 
formula predictions for the antialigned and nonspinning 
mergers. This may result from our failure to account for 
the energy carried away by gravitational-wave emission 
after the merger. While this emission becomes negligible 
compared to that during the inspiral as g — > 0, significant 
emission after the formation of a common horizon can 
occur for equal-mass mergers. We plan to incorporate 
radiation during the plunge and ringdown stages of the 
merger in future work. The final mass has not been reli- 
ably determined for the maximally spinning case, which 
perhaps accounts for the greater discrepancy between the 
different fitting formulas. These formulas included only 
a limited number of highly spinning mergers in the set 
with which they were calibrated. 

It is interesting to note when comparing rows K and 
KLP that superradiant scattering actually increases the 
predicted final mass for the maximally spinning merger, 
unlike in the other two cases. This is counterintuitive, 
since the top panel of Fig. [2] indicates that superradi- 



ant scattering extracts the most energy -BsR(Xjnsco) 
from the most highly spinning black holes. However, as 
the test particle inspirals, superradiant scattering reduces 
the spin x, which moves the ISCO radius rigco(x) out- 
wards and increases the ISCO energy -Eisco(x)- For large 
spins, the steep increase in i?isco(x) with spin more than 
compensates for the additional energy £^sr(X: ^sco) ex- 
tracted, thus increasing the final mass Mf. 

Until numerical relativists can accurately simulate 
mergers with high spins and small mass ratios in an ac- 
ceptable amount of time, we will need to rely on a com- 
bination of fitting formulas and test-particle extrapola- 
tions. These two methods complement each other, and 
any model that attempts to make predictions through- 
out the entire region < g < 1 should take advantage of 
both approaches. The fitting formulas are most accurate 
predicting the result of mergers close in parameter space 
to the simulations with which they were calibrated. Cur- 
rently this consists mostly of comparable-mass mergers 
with q > 0.1. The test-particle extrapolations like the 
one proposed in this paper are most reliable for g ^ 1. 
We can readily modify the AEI formula [31] for aligned 
mergers 

Xf = X + XV{s4X + S5V + to)+v{'^V3 + t2ii + t-sT]'^) , (17) 



with 



Xi 



9^X2 



(18) 



to incorporate the results of this paper. The first step in 
this modification was already taken in Rezzolla et al. [34j , 
where it was recognized that setting the coefficient of the 
term proportional to ry equal to 2-\/3 would reproduce the 
test-particle prediction for nonspinning black holes in the 
absence of superradiant scattering of gravitational waves. 
We propose that all terms linearly proportional to 77 can 
be replaced by our result for dxf/dq from Sec. 



IV 



XKLP = X + -J^{X)V + {h + S5xW + hif 



(19) 



This eliminates two of the coefficients from Eq. ( 17) and 



guarantees that the formula behaves properly in the test- 
particle limit. The highly accurate equal- mass nonspin- 
ning simulation fT? can still be used to determine ^2 a-s 
in Barausse and Rezzolla [31j if desired. If the formula 
no longer has enough degrees of freedom to reproduce 
existing numerical results, additional terms proportional 
to highers powers of 77 and x can be added. These new 
terms will not affect the test-particle behavior, as ry — >■ 
in this limit, but can help to fit comparable-mass mergers 
where rj < 0.25. There also remains additional freedom 
in our choice of %. Although the choice given in Eq. ( 18 ) 
possesses the desired symmetry and limiting behavior 



^0 ^ x->Xi 

00 ^ X ^ X2 , 



(20a) 
(20b) 
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this choice is not unique and alternatives may prove more 
suitable. An iterative approach for determining % as a 
function of {g,XijX2} as in BKL [53] and Kesden [5S] 
should be explored as well. 

This same approach can be used to improve fitting 
formulas for the final mass. Tichy and Marronetti [32] 
proposed the fitting formula 



M 

M 



f = 1 + 4(m° - 1)7? + IQm'^^Tfixi + X2) , (21) 



where m^ and m°'^ are fitting coefficients. In the test- 
particle limit 77 — > q ^ this reduces to 



Our Eq. (Hal in this limit predicts 






^E, 



ISCO 



(Xi) --BsR(xi,nsco)- 1, (23) 



correctly capturing the spin dependence of the binding 
energy Sisco and extracted energy i^g^ on the larger 
black hole's spin. We can use this result to replace the 
term linear in rj in Eq. (21 ) to yield 



M, 



KLP 



M 



1 



d_ 
drj 



M, 



KLP 



M 



{x)v + i^m^^ifx ■ (24) 



Notice that we have also replaced the sum of the spins 
in the third term with x, as the sum did not possess the 
required limiting behavior of Eq. ( 20 ) . Additional terms 
proportional to higher powers in 77 can be added to im- 
prove the agreement with comparable-mass simulations 
without affecting the test-particle behavior. 



VI. DISCUSSION 

The superradiant scattering of gravitational waves 
emitted during the inspiral sets a fundamental upper 
limit xiim = 0.9979 ± 0.0001 on the spin a black hole 
may attain by accreting test particles on quasicircular 
equatorial orbits. For this limiting spin, the energy and 
angular momentum advected when the test particle falls 
from the ISCO combines with that extracted by superra- 
diant scattering to leave the spin x = SjM"^ of the black 
hole unchanged. Black-hole spins greater than xiim will 
be reduced by mergers with test particles on quasicircu- 
lar equatorial orbits, even though the physical angular 
momentum S and mass M will both increase. 

Jacobson and Sotiriou [3S] recently argued that small 
but finite-mass particles on highly hyperbolic orbits could 
spin black holes above xiim and even above the Kerr limit 
itself. However, this study neglected gravitational radia- 
tion and any superradiant enhancement, which could be 
considerable for these highly spinning mergers. Such hy- 
perbolic orbits are also extremely unlikely in any realistic 
astrophysical scenario. 



Our limit is nearly indistinguishable from the limit 
= 0.9980 ± 0.0002 for black holes grown through 
thin-disk gas accretion [7] , suggesting that spin measure- 
ments alone cannot determine the relative contribution 
of these two channels of black-hole growth. If gas ac- 
cretion occurs through an advection-dominated accretion 
flow (ADAF) rather than a thin disk, the limiting black- 
hole spin is reduced to xadaf — 0.96, and the launching 
of jets can further reduce this spin to Xjots — 0.93 [TTj. 
Since our limit xiim is larger than for other modes of 
black-hole growth, it seems to be a fairly robust upper 
bound on spins that can be obtained by any astrophysical 
means. 

Although SBHs are grown primarily through gas ac- 
cretion rather than test-particle mergers, our limiting 
spin Xiim can still be astrophysically relevant. Gas ac- 
cretion occurs episodically, when galactic major merg- 
ers produce global torques that funnel large amounts of 
gas to galactic centers [35J [37]. Accretion of this gas 
onto the SBH leads to AGN feedback [35] which sim- 
ulations demonstrate is capable of suppressing further 
accretion and black- hole growth [5S]. By contrast, com- 
pact objects such as white dwarfs, neutron stars, and 
stellar-mass black holes (effectively test particles when 
orbiting an SBH) will inspiral continuously as they are 
scattered into the "loss cone" of orbits for which grav- 
itational inspiral occurs more rapidly than subsequent 
scattering [JO] . Although the same global torques driven 
by galactic major mergers may enhance the rate of these 
extreme-mass-ratio inspirals (EMRIs), other dynamical 
processes will refill the loss cone even between periods of 
AGN activity [Hj. Event rates for EMRIs depend sensi- 
tively on highly uncertain factors such as the dynamical 
state, star-formation history, binary fraction, and initial 
mass function within nuclear star clusters. The capture 
rate for lOM© black holes by a 10'' M© SBH can be as 
high as 10~ yr~^ immediately following a nuclear star- 
burst, and could average 10~^yr"^ for nucleated spiral 
galaxies like the Milky Way [12]. If such a galaxy has 
not had a major merger since redshift z ~ 1 (about 7.5 
Gyr ago), its SBH would aqcuire a fraction 



AM IOM0 X 10" Vr"^ X 7.5 Gyr 



Mn 



lO^M^ 



= 0.075 (25) 



of its mass through EMRIs. We see from Fig.[3]that this 
is a large enough fraction to drive the SBH spin towards 
our limiting value xiim from the initial spins xo <^ 0.9 
expected from gas accretion. This suggests that EMRIs 
can potentially affect SBH spin evolution, though modest 
gas accretion can also occur between major mergers and 
not all EMRIs will be on circular, equatorial orbits. 

Even if our limiting spin xum can be realized astrophys- 
ically, does its marginal difference from the Kerr limit 
have any observable consequences? The binding energy 
per unit mass at the ISCO 1 — -Eisco depends very sensi- 
tively on X near the Kerr limit, and sets an upper bound 
on the AGN efficiency e = L/Mc^, where L is the AGN 
luminosity, M is the accretion rate, and c is the speed of 



10 



light. Soltan [3D] argued that the total cosmological mass 
in SBHs could be estimated by equating the observed en- 
ergy emitted by AGN to that released during theoretical 
models of SBH growth. If spins at the Kerr limit were 
used for an analysis of this kind instead of our limiting 
spin xiim — 0.9979, the AGN efficiency e would be over- 
estimated by £;isco(xiim) - £^isco(l) = 0.1031 leading 
to a corresponding underestimate in the total SBH mass. 
In addition to their influence on the SBH spin, the EM- 
RIs themselves are an important source for gravitational- 
wave detectors such as LIGO and LISA. Decreasing the 
spin from the Kerr limit to xiim moves the ISCO radius 
out from risco(l) = rni to nsco(xiim) = 1.242toi, and 
decreases the ISCO orbital frequency from ^isco(l) = 
O.SOmj""^ to ^isco(xiim) = 0.42™^"^. This decrease could 
be observable in high signal-to-noise EMRIs seen by both 



these experiments. Our limiting spin xum is thus not only 
an interesting consequence of general relativity, but also 
one that has potentially observable implications for as- 
trophysics and gravitational-wave detection. 
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